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Abstract 

We exhibit a method for simultaneously treating recoil and threshold corrections in 
single-photon inclusive cross sections, working within the formalism of collinear factoriza- 
tion. This approach conserves both the energy and transverse momentum of resummed 
radiation. At moderate pr, we find the potential for substantial enhancements from higher- 
order perturbative and power-law nonperturbative corrections. 

Collinear factorization for short distance, inclusive cross sections is an important tool of 
particle physics. Next-to-leading order (NLO) corrections in the strong coupling a s are known 
for a wide, and still-growing, set of reactions. NLO cross sections, however, are not uniformly 
successful in describing hard-scattering data, which has suggested the necessity of examining 
yet higher orders in a s . Among the widely-discussed examples is prompt, or direct, photon 
production [p], Q]. 

In this letter, we describe a new approach for studying higher-order corrections in the context 
of collinear factorization for the prompt photon cross section, 

^ da AB ^ x (x T ) = ^ f dx ^ a/A ( Xa;;U ) f dx b <f> b / B {xb,v) p 3 j^zTZg for*/*) ^ (-q 
dp T ab J J ap T 
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where da a b^x / dpr is the hard-scattering function at fixed pr- Hadronic and partonic scaling 
variables are x T = 4:p T /S and x T = Ap T /s, respectively, with s = x a XbS the partonic center-of- 
mass energy squared. As s approaches its minimum value at x T — 1, the phase space available 
for gluon bremsstrahlung vanishes, which results in large corrections to do /dp? at all orders. 
Threshold resummation || f|, [| organizes this singular, but integrable, behavior of da/dpT- 

Another source of higher-order corrections is the recoil of the observed particle against unob- 
served radiation. Thus, for Eq. (|1|), beginning with NLO in da/dpr, a soft gluon may be radiated 
before the hard scattering, say a QCD-Compton process g q — > jq. The outgoing 7 q pair recoils 
against the soft gluon, and, as a result, the Compton process may be softer than would be the 
case without initial-state radiation. 

Valuable insights have emerged in studies of direct photon production in which the partons 
initiating the hard scattering are described by generalized parton distributions || [/J. Such dis- 
tributions include transverse momentum that is partly intrinsic and partly perturbative, whether 
formulated in impact parameter (b) space [Q or in transverse momentum space |], [10]. Never- 
theless, it is difficult to state confidently whether intrinsic transverse momentum is required by 



the data, which themselves allow varied interpretations [TT|. Li [113 showed first how to develop 
a joint resummation in both leading threshold and transverse momentum logarithms in parton 
distributions. 

Our approach remains within the formalism of collinear factorization. Contributions to the 
hard scattering function associated with threshold resummation are redistributed over soft gluon 
transverse momenta, simultaneously conserving energy and transverse momentum. Accounting 
for recoil leads to an additional enhancement. 

To compute higher-order recoil effects in Eq. ([!]), we consider the partonic cross section Q for 
ab — > 7 c, and fix the center-of-mass (cm.) rapidity, r\ of the photon. Near threshold, the overall 
process consists of a hard 2 — ► 2 subprocess, along with soft radiation, which can be factorized. 
To leading power, 1/Q T , this soft radiation does not change the flavor of the initial-state partons, 
and the hard-scattering subprocess recoils from soft radiation with a transverse momentum Q^. 
For fixed pr, the transverse momentum of the photon relative to the cm. of the hard scattering 
is p T — Qy/2. The cm. rapidity is thus related to Qr by 



4| P t-Qt/2| 2 _ 4|p^| 



cosh 2 77 



= x 2 T , (2) 



where s < S is the cm. energy squared of the 2 — > 2 hard subprocess. The relative transverse 
momentum, p^ thus sets the minimum value of s, and the scaling variable for the subprocess 
is x\. This kinematic linkage between transverse momentum and partonic energy drives the 
quantitative effects of recoil. We must limit the size of Qt in this analysis by some matching 
scale, to avoid going outside the region where the singularities in Qt dominate. For small Qt, 
however, we can jointly resum large corrections to da a b^ lc /d 2 QT dpT, in logarithms of 1 — x T 
and Qt- The latter cancel in the hard-scattering functions. Like threshold corrections, however, 
they may leave finite remainders. 

In summary, the resummed inclusive cross section for the leading contributions to prompt 
photon production is of the form: 



7 (resum)/— \ , (resum) 

Pt dp T ~ J d QT PT ePQ^dp-, 



, l resum ) 
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where p, is a matching scale. On the right-hand side, the direction of pt may be chosen arbitrarily, 
because of the azimuthal symmetry of the overall process. The inclusive partonic cross section 
Pr^ab^^c I ^Pt m Eq. ([3]) determines the resummed hard-scattering function p T da^^ /dpx in 
Eq. (P, after collinear factorization. 

The arguments leading to a jointly resummed formula for daj^^ /d 2 Q T dpr are similar to 
those for transverse momentum resummation in Drell-Yan production. The possibility of joint 
resummation for singular behavior in Qt and 1 — x\ is ensured by the factorization properties 



of the partonic cross section near threshold ||13|| , which we assume here. The resummed cross 
section that results from these considerations Hi is 



72o / C \ N+l 



ij " u \ 1 - X 



M e *- W (w) M^-lN' (4) 

in terms of moments of the physical, hadronic parton distributions and of the squared 2^2 
amplitude |M i:; | 2 , with defined in (0). The function Pij(N, Qt) is a "profile" of Qr-dependence 
for fixed N. Recoil is incorporated in Eq. (f|) through the interplay of the profile function and the 
factor (S/4p' T 2 ) N+1 . Singular Qt behavior is most easily organized in impact parameter space, 
where logarithms of both the transform variables b and N exponentiate. We may thus write the 
profile functions in Eq. (|) as 

Pa (N, Qt, Q,ft) = J d 2 b e- tb ' QT exp [E tJ ^ k (N, b, Q,n)\ , (5) 

where azimuthal symmetry insures that E is a function of b = |b| only. Here and below, 
Q = 2pt/xt denotes the hard scale in the exponent. The parameter fi represents both the 
renormalization scale and the factorization scale, whose explicit dependence we denote by 
below. Eq. (f|) reverts to a threshold-resummed prompt photon cross section when recoil is 
neglected, by setting b to zero in the exponents E of (||). 

We now construct the exponents Eij^k, to NLL in both b and N. As already noted by 
Li and Lai 0, logarithmic recoil corrections (logarithms of b) are generated entirely through 
initial-state radiation. Final state interactions, however, do produce logarithms of N. We may 
thus conveniently split each exponent into initial- and final-state parts: 

E^ k (N, b, Q, /i) = If (iV, b, Q, /i) + E™(N, Q, /i) . (6) 

In general, resummation of perturbative logarithms leads to nonperturbative power corrections 
in both b and N/Q. Because our primary interest is in transverse momentum distributions, 
however, we incorporate only a term proportional to b 2 in E ls ||, associated with the region of 
strong coupling [ 15 1 . The initial-state exponents E}^ will therefore be a sum of perturbative and 



nonperturbative contributions. The final state exponent, £" FS in Eq. (|6|) is essentially the same 
as for pure threshold resummation [|J. 

Rather than giving a formal derivation |TJ[ of the NLL perturbative exponent, we motivate our 



expression by requiring it to reproduce both threshold and kx resummations in the appropriate 
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limits. The following expression for the initial state exponent, including a nonperturbative term, 
gives all leading and next-to-leading logarithms in b and N: 



EfAN,b,Q,n) 



(2// 



A(a s (^ 2 )) + A> s (// 2 )) 



21n^ -b 2 F l3 (N,Q), (7) 



where N = Ne lE with 7^ the Euler constant. At this point, the coefficient Fij(N, Q) is arbitrary. 
The function A a (a s ) = Y. n ( a s/^) n A^ is, as usual fl6f , Aft = C a and 



Aft 
C A . 



-C a K = -C a 
2 2 



'67 _tt^ 
18 6 



10 



with C q = Cf, C g 

The perturbative exponent in Eq. (|7|) is similar to initial-state contributions in threshold 
resummation |5|], but with a new lower limit for the integral over the coupling scale //, which 
we denote Qx~ 1 {N,b). The following simple choice of x gives all NLL logarithms: 

x {N,b) = N + bQ/c 1 , (8) 

where C\ = 2e _TB . This doubly- resummed NLL exponent has the important property of respect- 
ing momentum conservation, not only for transverse components as in fc^-resummation, but for 
energy as well. 

The exponent in (|7]) is a function of the logarithmic variables, 



A 

P 



boQ> s (n 2 ) IniV 
6 a s (/i 2 )ln( X (iV,6)) 



(9) 



where b Q = (12ir) ^HC^ — 4T^Nf). The initial state contributions are given in these terms by 



Ej*(N,b,Q,»)= £ 



a=i,j L 

where, adopting the notation of Ref. || [10 

AW 



-^hft(\,(3)+hft(\,P,Q,Li,Li f ) 



(10) 



hft(X,P) 

hft(\,p,Q,v,fi f ) 
1 



2Kb 2 

Afth 
2ixbl 



[2/3 + (1 - 2A) ln(l - 2/3)] 



(11) 



-ln 2 (l-2/3) 



1-2A 

1-2/3 



(2/3 + ln(l-2/3)) 



27rfe n 



1 — 2A 

2/5 r^ + ln(1 ~ m 



A (i) 

-f-Aln 
tt6 



(12) 



with &! = (17C| - 10C A T R N F - 6C f T r N f )/24tt 2 . 

Given the above, we are in a position to illustrate the role of recoil, using Eq. (H), through the 
evaluation of the N and Q F integrals, with profile functions at fixed N defined by Eqs. (|5])-(12 



As in threshold resummation, it is necessary to define the integrals in a manner that avoids 
the ("Landau") singularities in the running coupling of Eq. (0) at large N and/or b. For pure 
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threshold resummation, a number of methods have been employed to define the N integral in 
djj) . We will use a "minimal" approach || [l7j , defining the contour C in (§) to extend to infinity 



through the negative real half-plane. For the b integral in the profile function, Eq. ([|), we have 
developed a new integration technique. 

To perform the b integral, we rewrite Eq. (|^) as 

poo 

Pij (N, Q T ,Q,/i) = 7r / dbb[h 1 (bQ T ,v)+h 2 (bQ T ,v)]e E ^^ N ' b ' Q ^. (13) 

J 

Here we introduce two auxiliary functions h 12 {z,v), related to Hankel functions and defined in 
terms of an arbitrary real, positive parameter v by integrals in the complex #-plane: 

hfav) = / d6e- tzsind 



h 2 (z,v) = / d6e~ lzsmy . (14) 

The hi t 2 become the usual Hankel functions Hi^(z) in the limit v — > oo. They are finite for any 
finite values of z and v. Their sum is always: hi(z, v) + h 2 {z^ v) = 2Jq(z), independent of v. The 
utility of the /i-functions is that they distinguish positive and negative phases in Eq. fllTf ), making 
it possible to treat the b integral of the profile function as the sum of two contours, one for each 
hi. These contours avoid the Landau pole by a deformation into either the upper half-plane (hi), 
or the lower half-plane (^2)- Such a definition of the integral is completely equivalent to the 
original form, Eq. (|5]) when the exponent is evaluated to finite order in perturbation theory. It 
defines the resummed integral "minimally" ||17| , without an explicit cutoff. 

Proton-nucleon cross sections computed in this fashion are illustrated in Fig. |I[ Here we 
show, for several values of photon p?, da^^x/ 'dQrdpr, the distribution of the cross section in 
the recoil momentum Qt- The kinematics are those of the E706 experiment 0. Since this is a 
"demonstration" calculation, we pick a nominal value of Fij(N,Q) = 0.5 GeV 2 for the Gaussian 
coefficient in Eq. ([7]), independent of parton type. We leave for future work a more realistic 
determination of this coefficient, including its Q-dependence. The parton distributions are those 



of Ref. [I8|, and we treat NLO iV-independent ("hard virtual") terms as in pi. For simplicity, 



we approximate Q in Eq. by 2px- Finally, we set v — 1 in Eq. ([14]) . 

The dashed lines are dap^^x / dQrdpr from Eq. @, but with recoil neglected by fixing 
S/4p' T 2 at S/4pt- The dashed lines thus show how each Qt contributes to threshold enhance- 
ment. The profiles are subtracted in their 6-dependent one- loop corrections, which produce the 
peaks and dips at low Qt for the smaller values of pt- Every curve shows a peak (near 2 GeV) 
associated with resummation, and falls off as a power with increasing Qt- Integrating the dashed 
lines reproduces the threshold-resummed cross sections of Ref. ||, after a subtraction at order 
a s to recover exact NLO. 

The solid lines of Fig. |l] show the iV-integrated distributions in Qt-, 
d<JpN™x / dQrdpr, now found by including the true recoil factor (S/Ap' T 2 ) N+1 in Eq. ([J). These 
curves therefore describe the fully resummed cross section. The resulting enhancement is clearly 
substantial. For small pr, the enhancement simply grows with Qt (it must diverge at Qt = %Pt), 
while for pt above 5 GeV it has a dip at about Qt = 5 GeV, which becomes more pronounced as 
Pt increases, but which never really gets to zero. This behavior makes it problematic to identify 
a unique matching scale /1 in (HI) that gives a stable prediction for the cross section. 
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This reservation notwithstanding, and reemphasizing that our calculation is primarily an il- 
lustration, rather than a quantitative prediction, we plot in Fig. |2] the resummed cross section 
for pt > 3.5 GeV, with the choice p, = 5 GeV for the calculation of Fig. [I], using the approximate 
procedure of Ref . || to convert do- pN ^ jX / dpr to a cross section integrated over a finite rapidity 
interval. We include an NLO photon fragmentation component in the cross section, again cal- 
culated as in ||]. To show the size of the enhancement that recoil can produce, as well as its 
potential phenomenological impact, we also exhibit the pure threshold resummed cross section 
and the E706 direct photon data in this range f2|. For the theoretical curves, we choose \x = pr', 
both resummed curves have sharply reduced factorization scale dependence compared to NLO. 
Evidently, recoil can be phenomenologically relevant. 

It will take some time to explore this resummation formalism, including the implementation 
of practical nonperturbative estimates and of matching procedures. Nevertheless, we hope that 
this formulation of recoil effects at higher orders, in the language of collinear factorization, is a 
step toward clarifying what has been a thorny issue in the application of perturbative QCD. We 
also anticipate that the joining of threshold and transverse momentum resummation [|l^], along 
with the minimal evaluation of the profile function described above, will have useful applications 
to Drell-Yan and related cross sections. 
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Figure 1: The prompt photon cross section da p ^^x / dQrdpx at \fs = 31.5 GeV, as a function of 
Q T for various values of photon p T . Dashed lines are computed without recoil (p^ = pr in ([|)), 
solid lines are with recoil. 
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Figure 2: Prompt photon cross section Ed 3 <T pN ^ lX /dp 3 for pN collisions at ^Js = 31.5 GeV. 
The dotted line represents the full NLO calculation, while the dashed and solid lines respectively 
incorporate pure threshold resummation |J and the joint resummation described in this paper. 
Data have been taken from |2j . 
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